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We show that finite angular momentum pairing chiral superconductors on the triangular lattice 
have point zeroes in the complex gap function. A topological quantum phase transition takes place 
through a nodal superconducting state at a specific carrier density x c where the normal state Fermi 
surface crosses the isolated zeros. For spin singlet pairing, we show that the second nearest neighbor 
d + id- wave pairing can be the dominant pairing channel. The gapless critical state at x c ~ 0.25 has 
six Dirac points and is topologically nontrivial with a T 3 spin relaxation rate below T c . This picture 
provides a possible explanation for the unconventional superconducting state of Na^CoC^ ■ yH^O. 
Analyzing a pairing model with strong correlation using the Gutzwiller projection and symmetry 
arguments, we study these topological phases and phase transitions as a function of Na doping. 
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The sodium cobaltate Na^OoC^, a layered triangular 
lattice electron system, has received widespread inter- 
est since the recent discovery of the 5K superconducting 
(SC) phase in water intercalated NazCo02 • 2/H2O near 
x = 0.3 [1]. While a very rich phase diagram is emerging 
for a broad range of Na concentrations [2| , the nature of 
the superconducting phase has remained poorly under- 
stood. The main evidence that the SC phase is uncon- 
ventional comes from the absence of the coherence peak 
in the in-plane NMR spin-lattice relaxation rate at T c 
and its power law temperature dependence (T 3 ) below 
T c 0, IHlH • Despite this evidence for an anisotropic SC 
gap function with line nodes, consistent with the spe- 
cific heat [f| and fxSR measurements, there has been 
considerable debate over the pairing symmetry. Measure- 
ments of angle-averaged Knight shift in powdered sam- 
ples have produced inconsistent results supporting either 
spin-triplet Q or spin-singlet pairing [j| . The most recent 
measurements on high quality single crystals [To] show 
that the spin contributions to the Knight shift decreases 
below T c along both the a and c-axis, which strongly sup- 
ports that the Cooper pairs are formed in the spin-singlet 
state. The Knight shift and the relaxation rate above T c 
show antiferromagnetic (AF) correlations (EL Hcj|. 

On the theoretical side, there has been a growing num- 
ber of proposals for unconventional superconductivity in 
the cobaltates. For spin-singlet pairing in a non-s-wave 
channel, the six-fold symmetry of the triangular lattice 
requires the low angular momentum paired state to have 
the chiral d x i_ y i ± id xy symmetry, raising the exciting 
possibility of a time-reversal symmetry breaking super- 
conductor. Earlier studies that drew analogy to the high- 
T c cuprates indeed found unanimously spin-singlet, d+id 
pairing via the nearest neighbor AF superexc hange in the 
triangular lattice t-J model at low doping (TlLll2Lll3,[T3|. 
However, it is conventional wisdom that a d + id paired 
state, as well as other chiral SC states with complex or- 



der parameters, has a full j?a_p and is thus inconsistent 
with NMR experiments [H, 0, HJ . 

In this paper, we show that extended chiral paired 
states beyond the nearest neighbor (NN) have generic 
point zeroes in the complex gap function inside the first 
Brillouin zone. Consider (£, n)-wave pairing with angu- 
lar momentum I for the pairs on the n-th NN bond, the 
order parameter can be written down in real space, 

A y = A fa e^, (1) 

where £ labels the angular momentum of the pair and 9ij 
is the angle of fij = f\ — fj between sites i and j. The 
Fourier transform of Ay is given by Ai n (k) — 2[/3L(fc) + 
i(3'L (fe)], where the real (/?') and the imaginary (f3") parts 
of the gap function are given in Table I in terms of the 
triangular lattice harmonics for the chiral p-wave (£ = 
1) and ci-wave (£ = 2) pairing. Nodes in the complex 
gap function arise where the real and imaginary parts of 
Afe(fc) vanish simultaneously, i.e. at the crossing points 
of the lines of zeroes of @'(k) and (3"(k). Fig. 1 shows 
the locations of the gap nodes in the first zone for chiral 
p and d-wave pairing. For first NN pairing, the nodes are 
pinned to the zone center and zone boundary and thus 
a generic chiral (£, l)-wave superconducting state has a 
full gap. Remarkably, for the case of n > 1, new nodes in 
the gap function appear inside the zone and can support 
nodal {I, n > l)-wave superconducting states. 

Hereafter we focus on the spin-singlet d+id- wave case, 
i.e. chiral (2,n)-paired states and argue that the exis- 
tence of the nodes is relevant for understanding the super- 
conducting state in hydrated cobaltates. For the second 
NN pairing with n — 2, there are six nodes inside the first 
zone marked by the solid circles in Fig. lb. They are lo- 
cated at ±k*, with a = 1,2,3, and k| = 2tt[2/(3%/3), 0], 
= 2tt[1/(3\/3),1/3], and k^ = 2tt [1/(3^3), -1/3], co- 
inciding with the six corners of the hexagonal y/3 x \/3 
reduced zone boundary. 
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Chiral p-wave pairing (£ = 1): 

n fl,„(fc) /%„(fc) 

1 \/3 sin ^ k x cos | k y sin fc M + cos ^k x sin 

2 — y^sin cos -^fco; sin ^3/^ + cos sin -^fc^ 

3 \/3 sin yJHk x cos fc a sin 2fcjy H~ cos sm Ky 



Chiral d-wave pairing (£ = 2): 

n /%,„(fc) /%„(*=) 

1 cos fc a — cos -^fea; cos ifcy %/3 sin fc^ sin i k y 

2 cos v^fc^ — cos §fc H cos -^fcx — \/3 sin sin -^fca; 

3 cos 2k y — cos \/3fea: cos ky \/3 sin \/3&2; sin k y 

TABLE I: Gap function A tn (k) = 2 [$„(fc) + if3"„(k)]. 




l=\ • 2nd-NN x = 0.25 hex. FS / = 2 

FIG. 1: Nodes of chiral [l, n)-wave pairing where the gap 
function A^ n (fc) vanishes, (a) p-wave and (b) d-wave with 
normal state FS at x — 0.25 (dashed line). 



The first indication that this is special for the trian- 
gular lattice cobaltates conies from the fact that a single 
hexagonal Fermi surface (FS) at the SC concentration 
x = 1/3 matches the V3 x %/3 reduced zone boundary 
such that these nodes would lie directly on the FS. Re- 
cent angle-resolved photoemission (ARPES) experiments 
in the both the unhydrated [HI, El an d the hydrated 
17| compounds near x ~ 0.3 indeed observe a single 
rounded hexagonal FS that coincides well with the hexag- 
onal %/3 x y/3 reduced zone boundary. This is highly 
unexpected since the Co 4+ has 5 (i-electrons occupying 
the three t 2g orbitals. It turns out that strong corre- 
lation effects beyond the band theory renormalize the 
crystal field splitting and bandwidths, leading to a single 
quasiparticle band crossing the Fermi level as observed 
by ARPES experiments [18|. In Fig. lb, we show the FS 
calculated in Ref. [l8| at x c = 0.25, which passes through 
the six gap nodes for the second NN d + id- wave pairing. 
The value of x c is smaller than 1 /3 due to the rounding 
of the hexagonal FS in both the ARPES data [l5| and 
the theoretical calculations [l8(. Interestingly, third NN 
d + id-wave pairing introduces six nodes (Fig. lb) at the 
corners of the hexagonal 2x2 reduced zone, which in- 
tersects the FS at x = 1/2 where the cobaltate is in an 
insulating phase with charge and spin order 0, [l^ . 

The second, direct evidence for nodal d + id-wave pair- 
ing at x — x c comes from the most recent NMR ex- 
periments. In high quality single crystals, Zheng et. 
ol. discovered that the T 3 decay of the relaxation rate 
1/Ti, i.e. the signature of line nodes in the gap func- 
tion, is obeyed down to the lowest temperatures only at 
x c ~ 0.26, whereas deviations from T 3 are observed on 
both sides of x c Q . The existence of the Dirac nodes in 
the gap function at the critical doping x c for spin-singlet 
pairing on the triangular lattice strongly favors the sce- 
nario of second NN d + id-wave pairing. In the following, 
we report a variational study of the (2, n)-wave pairing 
states in an effective single-band t — U plus pairing model 
of the cobaltates, calculate the Knight shift and the spin 
relaxation rate 1/Ti, and describe the topological prop- 



erties of the superconducting phases. 

The Hamiltonian is given on a triangular lattice, 



H 



hi 



where c\ a creates a hole of spin a and U is the Co on- 
site Coulomb repulsion. To describe the a\ g band cross- 
ing the Fermi level, we fix the first three NN hopping 
Uj = (h,t 2 ,t 3 ) = (-202,35,29) meV [19] . The hole den- 
sity rii — 1 — Xi where Xi is the electron doping concen- 
tration. The last term in Eq. @ is a phenomenological 
singlet pairing interaction Wij > between n-th NN sites 
i and j with Ay = CjfCjj — CjjCjj. The effects of strong 
correlation are accounted for in the Gutzwiller approxi- 
mation (20j , leading to a renormalized band dispersion, 

£fc = 2g t ti(cosk y + 2 cos \[Zk x j2 cos k y /2) 
+ 2g t t 2 (cos V3k x + 2 cos V3k x / 2 cos 3k v / 2) 
+ 2g t t 3 (cos 2k y + 2 cos V3k x cos k y ) — [i, (3) 

where gt = 2a;/ (1 + x) is the Gutzwiller renormalization 
factor in the large-U limit 0, H(J. The BCS state for 
(2, n)-wave pairing can be written equivalently as [24| 



i*n> = n 



u n k\ 1/2 exp ( -g n kc\^_ ki 



|0>, (4) 



where g nk = v nk /u nk = -{E nk - &, )/A* fc and the quasi- 
particle excitation energy E nk = + |A„fc| 2 . The gap 
function is A„ fe = 2W n A n [f3' 2n (k) + if3' 2 \ n {k)}, with /3 
given in Table I. We first study the variational ground 
state obtained by minimizing the energy with respect to 
the only variational parameter A„ for a given W n . In 
Fig. 2a, we compare A n for the 1st, 2nd, and 3rd NN 
pairing with identical pairing strength lFi.2,3 = 50meV 
at several dopings. The 2nd NN d + id pairing is the 
strongest with the largest pairing order parameter. In- 
cluding a NN Coulomb repulsion V 21| would further 
suppress the 1st NN pairing attraction W\, making the 
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FIG. 2: (a) Comparison of d + id pairing order parameters at 
several dopings for Wi,2,3 = 50 meV. (b-d) Properties of 2nd 
NN d + id- wave pairing at x — 0.19 = x c for Wi = 50meV 
and x = 0.4 > x c for W 2 = 73meV: Tunneling DOS (b), 
temperature dependence of Knight shift K s (c) and NMR 
relaxation rate 1/Ti (d) normalized by their values at T c . 



2nd NN d + id-wave the dominant spin-singlet pairing 
channel that drives the SC transition. We emphasize 
that while the subdominant channels can emerge below 
the transition, the admixture of a small component with 
n^2 would only result in small shifts of the six- fold sym- 
metric point nodes in the gap function and not change 
the results discussed here. 

Next we turn to the properties of the second NN d + id 
state. For our band parameters, the FS is almost cir- 
cular and x c ~ 0.19. At x = x c , the normal state FS 
crosses the six nodes at ±k* around which the quasi- 
particle dispersion Ek (hereafter we drop the index n) 
has a conical spectrum. For example, expanding around 

±k;, E k ~ yJ(A* + Bi)[k x T 4tt/(3V3)] 2 + S2fc|, with 
A = 3(ti - 2t 3 ) and B = 9WA/2. Thus, the six nodes 
are described by three pairs of Dirac fermion doublets 
with anisotropic velocities. They govern the properties 
of the low energy excitations. We calculate the tunneling 
density of states (DOS) N(E) and the temperature de- 
pendence of the a6-plane Knight shift (K s ) and the spin 
relaxation rate (1/Ti) according to 



K s , y x > x - 



{l,TN(E)}N(E)^jp-dE. (5) 



The results are plotted in Figs. 2b-d. The calculated 
N(E) shows a linearly vanishing, V-shaped DOS at x = 
x c generic of nodal d-wave pairing. As a consequence, 
the Knight K s and 1/Ti follow the T and T 3 behaviors 
below T c in agreement with the NMR experiments on 
single crystals at x ~ 0.26 [!, [l(|. For x ^ x c , the 
normal state FS does not overlap the nodes in the gap 
function. Fig. 2b shows that at x = 0.4 the DOS is fully 



gapped at low energies, but reverts back to the V-shape 
above the energy gap. This leads to a rapid crossover of 
K s and 1/Ti from the power law behaviors just below 
T c to exponential decays at low temperatures (Figs. 2c 
and 2d). In this case, we expect that disorder- induced 
filling of the DOS gap to produce crossovers to constant 
K s and linearly vanishing 1 /T\ at low temperatures Q . 

We now discuss the topological properties of the chi- 
ral SC states and the topological quantum phase transi- 
tion as the carrier density x moves across x c . The topo- 
logical order in chiral fermion-paired states has been a 
subject of growing interests in connection to unconven- 
tional superconductors /superfluids and fractional quan- 
tum Hall states 
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The key point is that the 
complex order parameter A/, and the dispersion form 
a unit pseudospin vector introduced by Anderson [25| : 
m(/c) = (ReAfe, — ImAfe, ^)/i?fc. Since |m| 2 = 1, m lives 
on a 2-sphere S 2 . The SC state can be viewed as a BCS 
mapping from the k-space, which is a torus T 2 (compact- 
ihed from an infinite plane) to the pseudospin S . This 
becomes evident when we write m in term of the pairing 
function g k = Vk/uk in Eq. Q, 



mi(fc) + irri2{k) 



> TO s(fc) = T~r i — fa ■ ( 6 ) 



i + k-l 



i + kl s 



Such maps are classified by homotopy classes with topo- 
logical winding numbers Q. The lattice version of Q is, 



*2 = g~ m ( fcl ) ' [ m ( fc 2) x m(fc 3 

A 



(7) 



where the sum is over all elemental triangular plaque- 
ttes with corners labeled as 1, 2, 3. Topologically distinct 
phases are categorized by different winding numbers Q 
that count the number of times the pseudospin configu- 
ration m wraps around the sphere. We find that Q ^ 
if gk has a singular vortex structure near the nodes of 
the complex gap function, a condition that depends on 
the location of the FS. In Fig. 3, we show the calculated 
Q for x > x c , x < x c , and at the critical state x = x c . 
All of them are topologically nontrivial. For x > x c , 
the normal state FS encloses only the gap node at the 
zone center where gk oc l/(k x —ik y ) 2 , exhibiting a singu- 
lar double vortex. This charge-two magnetic monopole 
contributes two flux quanta to the surface integral and 
gives Q = q=2 (Fig. 3) for the fully gapped topological 
phase. The sign corresponds to the relative sign of the 
real and imaginary parts of the order parameter. Note 
that near the nodes outside the FS, gk behaves as zeroes, 
e.g. gk oc k x =p 47r/ (3\/3) + ik y near ±k*, and thus do not 
contribute to Q. As the doping is reduced to x < x c , the 
normal state FS encloses, in addition to the node at T, 
the six nodes at ±k* where gk exhibits singular vortex 
structure; e.g. gk oc l/[k x ^ 47r/(3vo) + ik y ] near ±k*. 
Each pair of the nodes can be viewed as a singular double 
vortex. The contribution to the magnetic flux through 
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FIG. 3: Topological winder number as a function of doping 
in the second NN d + id-wave superconducting state. Insets 
show the locations of the nodes in the complex gap function 
and the normal state FS. 



the sphere is thus increased by six flux quanta leading to 
Q = ±4 as shown in Fig. 3. 

In the critical state at x = x c , the FS passes through 
the six gap nodes. The pseudospin vector m is not well 
defined exactly at these isolated diabolic points, which 
are lattice- regulated in the summation in Eq. ([7]). Re- 
markably, Fig. 3 shows that the winding number Q = ±1 
and the critical state is a gapless topologically nontrivial 
state. This suggests that each pair of the gap nodes in 
the critical state contributes a single vortex in fc-space , 
reminiscent of the situation in the p 
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ip state 

Indeed, evaluating Q in the continuum limit along the 
line integrals of the six circles cutting out the nodes at 
k* , one finds that each of the latter contributes half a 
magnetic flux leading to the total Q = ±1. We note that 
the existence of such a gapless topological phase in two 
dimensions is new and consistent with the view point that 
gapless fermions emerge as a result of a form of quantum 
order [27j in the critical state. Physically, the topolog- 
ically invariant winding number Q corresponds to the 
quantization of the spin Hall conductance a xy [H, [24| • 
In unit of the spin conductance quantum (h/2) 2 /2Trh, 
the spin Hall conductance is a xy = Q. Thus, we predict 
that as the electron doping x evolves across x c , the super- 
conducting state of the cobaltate exhibits quantum spin 
Hall transitions from =p2 to ±4 with quantized critical 
spin Hall conductance ±1. We expect novel properties 
associated with the edge states. 

In summary, we have proposed a class of nodal chiral 
superconductors. The spin-singlet, second NN d + id- 
wave pairing turns out to be consistent with the NMR 
experiments on high quality single crystals of hydrated 
cobaltate superconductors, and has some remarkable 
topological properties related to the quantum spin Hall 
transitions. Although the microscopic origin for the sec- 
ond NN d + id pairing is beyond our scope here, it is 
conceivable that (i) a strong NN Coulomb repulsion [2l| , 
(ii) a larger 2nd NN superexchange, and (iii) the frus- 
tration of the NN antiferromagnetic correlation on the 



triangular lattice and the proximity to inhomogeneous 
charge/spin ordered state Qjl, Hl| will favor a dominant 
extended pairing interaction beyond the first NN. It is 
noted that the chiral pairing state breaks time-reversal 
symmetry which should in principle be detectable by /iSR 
or optical Kerr experiments with sufficient resolution. 
The orbital current near a unitary impurity produces a 
static magnetic field. To estimate the size of the field, 
we performed self-consistent calculations for the circu- 
lating current near a unitary impurity. The maximum 
current around the impurity site is ~ lOOnA for a pairing 
strength W% — 50meV corresponding to a T c ~ 5K . This 
gives an estimate of the magnetic field at the impurity 
site ~ 1G, which is close to the upper bound set by the 
earlier /xSR experiments 28, 2^]. Future experiments are 
very desirable to determine whether time-reversal sym- 
metry is broken in the superconducting cobaltates. 

We are grateful to Y. Yu for many useful discus- 
sions. This work was supported by DOE grant DE-FG02- 
99ER45747 and NSF grant DMR-0704545. 



[1 
[2 
[3 
[4. 
[5. 

[r\ 

[8 

[9 
[10 

[11 
[12 

[13 
[14 

[is; 

[16 

li- 
tis 

[19 
[20 

[21 

[22 
[23 

[24 
[25 
[26 
[27 
[28 
[29 



K. Takada, et al, Nature (London) 422, 53 (2003). 
M. L. Foo, et al, Phys. Rev. Lett. 92, 247001 (2004). 
T. Fujimoto, et a/.,Phys. Rev. Lett. 92, 047004 (2004). 
K. Ishida, et al, J. Phys. Soc. Jpn, 72, 3041 (2003). 

G. -q. Zheng, et al, J. Phys.: Condens. Matter 18, L63 
(2006). 

H. D. Yang, et al, Phys. Rev. B 71, 020504 (2005). 

A. Kanigel, et al, Phys. Rev. Lett. 92, 257007 (2004). 
M. Kato, et al, J. Phys.: Condens. Matter 18, 669 
(2006). 

Y. Kobayashi, et al, J. Phys. Soc. Jpn. 74, 1800 (2005). 
G.-q. Zheng, et al, Phys. Rev. B 73, 180503 (2006). 

G. Baskaran, Phys. Rev. Lett. 91, 097003 (2003). 

B. Kumar and B. S. Shastry, Phys. Rev. B 68, 104508 

(2003) . 

M. Ogata, J. Phys. Soc. Jpn. 72, 1839 (2003). 

Q.-H Wang, D.-H Lee, and P. A. Lee, Phys. Rev. B 69, 

092504 (2004). 

H. -B. Yang, et al, Phys. Rev. Lett. 95, 146401 (2005); 
Phys. Rev. Lett. 92, 246403 (2004). 

M. Z. Hasan, et al, Phys. Rev. Lett. 92, 246402 (2004). 

T. Shimojima, et al, Phys. Rev. Lett. 97, 267003 (2006). 

S. Zhou, et al, Phys. Rev. Lett. 94, 206401 (2005). 

S. Zhou and Z. Wang, Phys. Rev. Lett. 98, 226402 (2007). 

D. Vallhardt, Rev. Mod. Phys. 56, 99 (1984); F.C. Zhang, 

et. al. Supercond. Sci. Technol. 1, 36 (1988). 

O. I. Motrunich and P. A. Lee, Phys. Rev. B 69, 214516 

(2004) ; 70, 024514 (2004). 

G.E. Volovik, JETP Lett. 66, 522 (1997). 

T. Senthil, J.B. Marston, and M.P.A. Fisher, Phys. Rev. 

B 60, 4245 (1999). 

N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
PW. Anderson, Phys. Rev. 110, 827; 112, 1900 (1958). 
G.E. Volovik, Sov. Phys. JETP 67, 1804 (1988). 
X. G. Wen and A. Zee, Phys. Rev. B 66, 235110 (2002). 
W. Higemoto, et al, Phys. Rev. B 7 0, 134508 (2004). 
Y. J. Uemura, et al., |arXiv:cond-mat/0403031| (2004). 



